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ABSTRACT 

We propose a new approach to the missing baryons problem. Building on the common assumption that the missing baryons are in 
the form of the Warm Hot Intergalactic Medium (WHIM), we further assumed here that the galaxy luminosity density can be used as 
a tracer of the WHIM. The latter assumption is supported by our finding of a significant correlation between the WHIM density and 
the galaxy luminosity density in the hydrodynamical simulations of Cui et al. (2012). We further found that the fraction of the gas 
mass in the WHIM phase is substantially (by a factor of — 1.6) higher within the large scale galactic filaments, i.e. ~70%, compared 
to the average in the full simulation volume of ~0.1 Gpc 3 . The relation between the WHIM overdensity and the galaxy luminosity 

overdensity within the galactic filaments is consistent with linear: <5 w him = 0.7 ± 0.1 x <5^ ±0 ' 2 - 

We applied our procedure to the line of sight to the blazar H2356-309 and found evidence for the WHIM in correspondence of the 
Sculptor Wall (z ~0.03 and log N H = 19.9/(q|) and Pisces-Cetus superclusters (z ~0.06 and log Nh = 19.7/g 2 ), in agreement with the 
redshifts and column densities of the X-ray absorbers identified and studied by Fang et al. (2010) and Zappacosta et al. (2010). This 
agreement indicates that the galaxy luminosity density and galactic filaments are reliable signposts for the WHIM and that our method 
is robust in estimating the WHIM density. The signal that we detected cannot originate from the halos of the nearby galaxies since 
they cannot account for the large WHIM column densities that our method and X-ray analysis consistently find in the Sculptor Wall 
and Pisces-Cetus superclusters. 

Key words. Cosmology: observations - large-scale structure of Universe - intergalactic medium 


1. Introduction 

A large fraction of the local (z < 1) baryons as predicted by the 
concordance ACDM cosmology are not detected (e.g. Nicastro 
et al. 2008; Shull et al. 2012), giving rise to the missing baryons 
problem. Since the current cosmological model agrees very ac¬ 
curately with most of the observations, it is likely that the prob¬ 
lem of these missing baryons is observational by nature, rather 
than a problem with the cosmological theory. 

In the current view the cosmological large-scale structure 
formation is driven by dark matter, due to its mass dominance 
and solely gravitating nature. First the dark matter skeleton of the 
cosmic web is formed, and then the baryons fall and condense to 
these filamentary gravitational potential wells (e.g. Cen & Os- 
triker 1999; Cui et al. 2012, hereafter C12). With cosmological 
time these dense filamentary regions become more pronounced 
and finally at z < 1 the released gravitational potential is suffi¬ 
cient to shock-heat the baryons at over-densities 5b = PbliPb ) in 
the range 1-100 to temperatures of 10 5 —10 7 K (e.g. Dave et al. 
2001; Bregman 2007). The surrounding medium is also enriched 
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by metals through galactic outflows and stellar feedback (e.g. 
Cen & Ostriker 2006; Schaye et al. 2010). The gas in this phase 
is usually called the Warm Hot Intergalactic Medium (WHIM). 

The expected low density of the WHIM renders it very dif¬ 
ficult to detect in emission and thus a viable candidate for the 
(so far) missing baryons. The rare detections of WHIM in emis¬ 
sion have been obtained by e.g. Kull & Bohringer (1999) and 
Zappacosta et al. (2005) using ROSAT/PSPC observations of the 
Shapley and Sculptor superclusters, respectively and by Werner 
et al. (2008) for a cluster pair A222/A223. The full understand¬ 
ing of these detections requires measuring the redshift of the 
emission, which is currently beyond the capabilities of Chandra 
and XMM-Newton telescopes. The situation will only improve 
with advent of Athena mission. The absorption measurements in 
X-rays (e.g. Zappacosta et al. 2010, hereafter Z10; Fang et al. 
2010, hereafter F10; Ren et al. 2014) and especially in far ultra¬ 
violet (Shull et al. 1998, 2003; Tilton et al. 2012; Stocke et al. 
2014) have been quite successful. However, they are limited by 
the small number of suitable background sources observable in 
the bright state behind significant WHIM structures. 
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The filamentary regions with enhanced gas density along the 
cosmic web, in addition to hosting WHIM, are also preferen¬ 
tial sites of galaxy formation and stellar feedback process, see 
(e.g. Baugh 2013; Cui et al. 2012; Cen & Ostriker 2006; Nusser 
et al. 2014; Sigad et al. 2000). Thus galaxies and the WHIM 
are expected to trace the same underlying large-scale filaments 
and, therefore, should be spatially correlated. Building on this 
assumption we have developed a method to trace the missing 
baryons in the form of WHIM using the observed galaxies and 
their luminosity density as a proxy. The method consists of 1) de¬ 
tecting the galaxy filaments in spectroscopic galaxy catalogues, 
2) determining the galaxy luminosity density fields around the 
detected filaments, 3) deriving and calibrating a relation between 
the galaxy luminosity density and the WHIM density from hy- 
drodynamical N-body simulations in cosmological boxes and 4) 
applying the above phenomenological relation to estimate the 
WHIM density using the observed luminosity density of galax¬ 
ies in the filaments. 

We tested the accuracy of our WHIM column density estima¬ 
tion by applying it to the 2dF data around the Sculptor Wall (S W) 
and Pisces-Cetus (PC) superclusters and comparing our values 
with those obtained via the X-ray absorption measurements for 
these structures (Buote et al. 2009; F10; Z10). 

Our definition of WHIM is the gas in the temperature 
range of 10-10 7 K and in the baryon overdensity range db = 
Pb/(Pb) of 1-100. We use O m = 0.3, Oa = 0.7 and H 0 = 
70 km s^Mpc -1 . 

2. Methods 

2.1. Filament detection 

The detection of filaments is performed by applying an object 
point process with interactions i.e. the Bisous process (Stoica 
et al. 2005) to the spatial three-dimensional distribution of galax¬ 
ies. The method provides a quantitative classification that agrees 
well with a visual impression of the filamentary nature of the 
cosmic web and is based on a robust and well-defined mathe¬ 
matical scheme. More details regarding the Bisous model can be 
found in Stoica et al. (2007, 2010) and Tempel et al. (2014a). For 
reader convenience, a brief summary is provided below. 

The Bisous model approximates the filamentary web by a 
random configuration of small cylindrical segments. The model 
assumes that locally galaxy distribution can be probed with rela¬ 
tively small cylinders, which can be combined to trace a galaxy 
filament if the neighboring cylinders are oriented similarly. One 
of the advantages of such approach is that it relies directly on 
the observed positions of galaxies and does not require any addi¬ 
tional smoothing kernels for creating a filamentary density field. 

The solution provided by the Bisous process is stochastic. 
It is thus expected that there is some variation in the detected 
filamentary patterns for different Markov-chain Monte Carlo 
(MCMC) runs of the model. However, thanks to the stochastic 
nature of the model, we gain a morphological and a statistical 
characterisation of the filamentary network simultaneously. 

In practice, we first fix the width scale of the filaments to 
1.4 Mpc since as found by Tempel et al. (2014a) this is the typ¬ 
ical scale of the filaments in the cosmic web as traced by galax¬ 
ies. The results are not very sensitive to the exact value of the 
scale. By applying the algorithm we then obtain a three dimen¬ 
sional “visit” map containing information related to the filament 
detection probability (see Tempel et al. 2014a, for the definition 
of this parameter). At the moment we do not have a method to 
relate this information to an exact detection probability. Rather, 


we set experimentally the lower limit for the visit parameter to 
0.005 since this choice removes the noisy outer regions of the 
filaments and since the higher values would break filamentary 
structures clearly identified by eye inspection into disconnected 
sub-units (see Section 3.3). 

2.2. Luminosity density field 

The goal of the luminosity density (LD) method is to build up a 
three dimensional luminosity density field using the spatial dis¬ 
tribution and the luminosities of a given galaxy sample (see Li- 
ivamagi et al. 2012, for details). In case of the observational data 
(in contrast to the simulated data) we start from the fluxes of 
galaxies which we convert into luminosities (including the K- 
corrections). 

Since the 2dF galaxy catalogue is flux limited (to bj ~ 19.5) 
the mean galaxy number density and luminosity decrease with 
the redshift. To account for this selection effect we assign a sta¬ 
tistical weight to each galaxy proportional to the fraction of the 
galaxy luminosity accessible to observations at a given redshift 
(see Tempel 2011; Tempel et al. 2014b). 

We then smooth the galaxy luminosity distribution (observed 
or simulated) with a symmetric B3-spline kernel function. The 
smoothing scale sets the minimum physical scale of the struc¬ 
tures that we can identify in our analysis. In this case we set it 
equal to 1.4 Mpc to match the characteristic scale of the fila¬ 
ments in the galaxy distribution (see Tempel et al. 2014a). 

We then sample the smoothed LD distribution at points of a 
uniform cubic grid with a grid size 1.4 Mpc 1 , encompassing the 
survey volume, thus creating the LD field. A smaller sampling 
scale would be desirable, but computationally expensive in case 
of large simulations (see below). Our adopted sampling scale of 
1.4 Mpc represents an acceptable compromise between the res¬ 
olution and the computational cost. Note that when we analyse 
the individual sight-lines in the observational data we increase 
the sampling resolution in order to derive smooth distribution of 
the luminosity density (see Section 4). 

3. Large-scale structure simulation 

As discussed in Introduction, numerical simulations of the build¬ 
up of the cosmic structures, as well as theoretical arguments, in¬ 
dicate that both galaxies and the WHIM trace the same underly¬ 
ing distribution of dark matter. The relation between the baryon 
tracers and the underlying mass is not deterministic since it also 
depends on the type of tracer, the scale, the underlying density 
and the ill-known stellar feedback processes. It also depends on 
time. In our analysis we consider a rather local sample, so the 
latter dependency can be ignored. 

However, in the filamentary regions, i.e. on scales signifi¬ 
cantly larger than those affected by stellar processes, the LD- 
WHIM density relation may be usefully tight. Building upon this 
assumption, we shall consider the hydrodynamical simulations 
carried out by C12 to derive a relation between the galaxy LD 
and the WHIM density within cosmic filaments. The procedure 
to build up and calibrate this relation is described in this Section. 

3.1. Gas and dark matter 

The simulations of C12 have been carried out using the TreePM- 
smoothed particle hydrodynamics (SPH) code GADGET-3, an 

1 In order to avoid under-sampling, the sampling scale should not ex¬ 
ceed our adopted smoothing scale. 
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improved version of GADGET-2 (Springel 2005). It follows the 
evolution of a box with a co-moving size of 590 Mpc per side 
considering both dark matter and gas. Each component is de¬ 
scribed using 1024 3 particles from z = 41 to the present epoch 
(see C12 for full details). 

We used these particle data as an input to our analysis. We 
divided the full simulation volume into a cubic grid whose grid 
size was determined from two criteria: the minimum gas over¬ 
density that we want to sample and the need to avoid oversam¬ 
pling the gas density field. Since we focus on the WHIM, we 
want to trace this gas down to a density contrast c>b = 1- This cor¬ 
responds to an inter-particle separation of 0.6 Mpc in the simula¬ 
tion, which is thus the minimal useful grid size. In order to avoid 
the oversampling of the field we set the grid size to 1.4 Mpc, 
for both the gas and the DM. With this choice, the gas and DM 
density field of the simulation are defined on a 410 3 cubic grid 
which constitutes the input data set to model the gas properties. 

The physical modelling of the gas component includes radia¬ 
tive cooling, star formation and feedback from supernova rem¬ 
nants, but not AGN feedback. Even if ignoring AGN feedback 
is rather unphysical, this assumption does not affect our work 
since the galaxy luminosity is self-consistently modelled in the 
framework of the halo occupation distribution model (see next 
Section) whereby galaxy properties are determined from the DM 
distribution on a statistical basis. In addition, since we shall fo¬ 
cus on a rather local galaxy sample within which evolutionary 
effects can safely be ignored, we only consider the z — 0 output 
of the simulation. 

We then created 3D maps of the gas density and tempera¬ 
ture and the dark matter density defined on a cubic grid as fol¬ 
lows. Each gas particle’s mass, m,-, is distributed over a variable 
number of cells with weight proportional to the SPH smoothing 
kernel integrated over the cell volume (see similar applications in 
Roncarelli et al. (2006, 2007, 2012). The gas density is then com¬ 
puted by dividing the mass associated to each cell by its volume. 
We then computed the mass-weighted gas temperature 3D-maps 
by smoothing the value of w,T, of each gas particle, where 7) 
is the SPH particle’s temperature, and then dividing each cell’s 
value by its mass obtained from the previous mapping. 

DM particles are treated differently. The DM density field is 
interpolated at the points of the same cubic grid as above, using 
the Clouds-In-Cells (CIC) method that assigns each DM particle 
to 8 cells with weight that depends on the particle position. 

3.2. Galaxies 

We created the galaxies for the above simulation (Cl 2) by popu¬ 
lating the DM halos using a halo occupation distribution (HOD) 
constrained with the SDSS data (Zehavi et al. 2011). In sum¬ 
mary, we used the masses and positions of each DM halo in or¬ 
der to obtain the magnitudes and positions of the galaxies as fol¬ 
lows (see Appendix A for more details): 1) in order to know the 
spatial scale of each halo, we defined its virial radius according 
to the spherical collapse formalism (see e.g. Peebles 1980; Eke 
et al. 1996; Kitayama & Suto 1996; Bryan & Norman 1998); 2) 
we populated each halo with subhalos, up to the virial radius, by 
performing a Monte Carlo realisation of the subhalo mass func¬ 
tion model (Giocoli et al. 2010); 3) for each halo we performed 
an abundance matching approach to assign a given central galaxy 
or a satellite with a luminosity using halo occupation distribution 
of Zehavi et al. (201 1). 

Our HOD implementation predicts galaxy luminosities in the 
r band, whereas we wish to apply the results to 2dF galaxies 
(see below) that were observed in the bj band. However, one 


Table 1. Conversion factors. 


name 

value 

<LD,.)“ 

0.0113 x 10 lu L o Mpc -3 

<LD,/ 

0.00772 x 10 10 L g Mpc" 3 

<Pb> c 

0.618 x 10 10 M q Mpc~ 3 

redshift CMB correction^ 

-9.4 x 10~ 4 


Notes. <a) The mean galaxy luminosity density of the full 2dF area in the 
bj band. (h) The mean galaxy luminosity density of the C12 simulation 
in the r band. (c) The cosmic baryon density at z = 0. (rf) To be added to 
heliocentric redshifts for conversion to CMB rest frame in direction of 
H2356-309. 


can extrapolate LD predictions from a band x to another band 
y , if the the mean LD has been estimated in both bands: as¬ 
suming that the luminosity overdensity (>ld = LD/(LD) for a 
given galaxy population is band-independent, one can estimate 
LD,, = LD t x (LD !/ )/(LD t ). The mean LD value at z = 0 in the 
/•-band obtained from the C12 simulation and that in the b j band 
obtained from 2dF (both estimated by us) are listed in Table 1. 

In addition, one should also consider the effect of morpho¬ 
logical segregation, i.e. the fact that red galaxies populate mostly 
the high density regions and blue galaxies are more homoge¬ 
neously distributed. Since we are interested in the WHIM fila¬ 
ments, i.e. we concentrate on the low to intermediate overdensity 
regions, our galaxy populations are expected to be dominated 
by blue galaxies. Indeed, based on our preliminary analysis, the 
morphological segregation of galaxies in filaments is very weak. 

3.3. Application of filament-finding and LD field methods to 
simulated data 

We constructed the r band luminosity density (LD,) fields by 
applying the LD method to the galaxy population created above. 
We first extracted the luminosity densities and gas densities in 
the full volume of the simulation studied in this work. We found 
that the fraction of the gas mass in the WHIM temperature and 
density range, i.e. the WHIM volume filling factor, is ~50%. 

We then applied the Bisous model to the galaxy distribution 
and defined volume elements with visit parameter values higher 
than 0.005 (see Section 2.1) as members of filaments. The ben¬ 
efit of this procedure is that it can be applied both to simulated 
and observational data, based only on the 3D distribution of the 
galaxies. However, in the case of the observational data, we do 
not have any temperature or density information of the stud¬ 
ied intergalactic gas. Thus, for consistency, we will not directly 
apply any temperature or density selection when extracting the 
simulated data. 

The environmental selection did in fact have the desired ef¬ 
fect of excluding the voids and the noisy, low temperature fila¬ 
ment edges and concentrating more on the filament cores (see 
Fig. 1). The WHIM mass fraction increased from 50% to 70% 
when we restricted the analysis to such environmentally-selected 
volume. In the following we thus assume that 70% of the gas 
mass in the regions satisfying the Bisous model filament detec¬ 
tion criteria is WHIM, i.e. p w hi m = 0.7 x p gas . 

We want to exclude galaxy clusters from the analysis since 
their high galaxy number density and low expected WHIM mass 
fraction (e.g. Branchini et al. 2009) yields a very different LD- 
WHIM density ratio from that in the filaments. The above fil¬ 
tering minimises rather efficiently the contribution of clusters of 
galaxies, proven by the fact that only ~10% of the selected gas 
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Fig. 1 . Gas density (10 10 M G Mpc 3 ) image of an interesting region from the simulation (02) when applying no filtering (left panel) and when 
applying the filamentary environment selection (right panel). The pixel size and the image thickness are both 1.4 Mpc. 


mass has temperature and densities that both exceed the WHIM 
upper limits. 

3.4. Determining the PDF of the WHIM density in different 
luminosity environments 

The relation between the WHIM density and the luminosity 
density, both stochastic tracers of the underlying mass den¬ 
sity field is, unsurprisingly, also stochastic. As a result, for a 
given LD, the values of the WHIM density spread over quite 
a broad range. Here we wanted to quantify more than the sim¬ 
ple spread: we characterised the constrained probability function 
■P(Pwhim I LD r ). 

We have restricted our analysis in the LD,. range of [0.01- 
2.0] x 10 l(, LoMpc~ 3 . This choice is justified on the basis of the 
simulations of C12 that indicated that ~ 95% of the WHIM mass 
of the full simulation volume is within this LD, range. Increasing 
the range would therefore decrease the “purity” of the sample. 
We divided this LD, range into 15 equally-spaced logarithmic 
bins and measured the conditional PDF in each of them. 

The frequency histograms of the WHIM gas density at a 
given LD, bin constitute our estimate of the conditional prob¬ 
ability distribution function P(f> xv i lml I LD,.) (see Fig. 2). The 
WHIM density value corresponding to the position of the maxi¬ 
mum correlates with the LD, value, which indicates the presence 
of a LD-WHIM density correlation that we explore in the next 
Section. As long as LD, is small the PDF is well approximated 
by a log-normal distribution. For Old > 10 a positive skewness 
appears that develops into a secondary peak at <5 w him ~ 30 which 
shifts to progressively larger p„hi m values when LD increases. 
These features suggest the presence of a population of luminous 
objects with associated WHIM gas with over-densities higher 
than ~ 10. We speculate that these objects correspond to lumi¬ 
nous galaxies or groups with relatively dense WHIM gas in their 
X-ray halos. The fact that the WHIM gas associated to this peak 
seems to follow the same LD-WHIM density relation as the gas 
in the main peak indicates that this is a real feature and not a 
simulation artifact. 

To test this hypothesis we went back to the simulations to 
characterise the mock galaxies and the gas that contribute to 


the secondary peaks. We found that the peak disappeared when 
we removed from the sample the data from volume elements 
with galaxies brighter than M r = -20.5 or gas colder than 
T = 5x 10 6 K (see Fig. 3). This indicates that the second peak 
is contributed by gas particles preferentially associated to bright 
galaxies in halos with masses in excess of a few times 10 12 M 0 
(see Appendix A). These objects are rarely isolated and instead 
are typically found in groups (see Appendix A). Thus, and con¬ 
sidering also the relatively low temperature, the secondary peak 
appears to be contributed by gas halos of galaxies within galaxy 
groups or by the intra-group gas, or both. 

3.5. LD - WHIM density relation 

Within the filamentary environments, the simulated WHIM den¬ 
sity and the LD, correlate well: the Pearson correlation coeffi¬ 
cient is ~0.80 while the number of data points is ~ 1.7 x 10 6 . 
This indicates that the LD traces well the WHIM. Thus, we 
utilise in the following this correlation for finding the WHIM. 

We obtained our main goal of estimating the relation be¬ 
tween LD and WHIM gas density as follows. 1) We considered 
the conditional probability functions P(p w i um | LD, ) in each one 
of the 15 LD,. bins. 2) We Monte Carlo sampled the individual 
PDFs to obtain a statistical realisation of the LD, - p w him rela¬ 
tion. 3) We fitted the corresponding data with a power-law rela¬ 
tion p W him = A x LDf and stored the best fit parameters A and B. 

4) We went back to step 1 and repeated the procedure 10 4 times. 

5) We estimated the resulting PDFs for A and B. 

For both parameters the PDF was well approximated by a 
log-normal distribution except in the positive tail in which the B 
parameters exhibit an excess probability (see Fig. 4). Using the 
centroids and the widths of the best-fit log-normal models to the 
distributions we thus obtained the LD, - p w hi m relation as 

Pwhim = 0.4 ± 0.1 x 10 10 M 0 Mpc- 3 x <5 ld °' 9±0 ' 2 (1) 

and an analogous one expressed in terms of gas overdensity 

d whi m — 0.7 + 0.1 x d LD °' 9±0 ' 2 , (2) 

where the luminosity overdensity <5 ld = LD/(LD) 

is the independent variable, <5 w him = Pwhim Kpb), and 
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-1.5 -1.0 -0.5 0.0 0.5 

log (Pwhim) [10'° M sun Mpc" 3 ] 





0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 
'og (Pwhim) [TO 10 M sun Mpc" 3 ] 




0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 
'og (Pwhim) [TO ’ 0 M sun Mpc" 3 ] 



0.8 1.0 1.2 1.4 1.6 1.8 2.0 
log (Pwhim) [’O ' 0 M sun Mpc" 3 ] 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 
log (Pwhim) [TO 10 M sun Mpc" 3 ] 


Fig. 2. The WHIM density distributions (crosses) and the best-fit log-normal models (lines) within selected LD bins. h LD indicates the central 
luminosity overdensity value of a given bin. (J whim indicates the WHIM overdensity value corresponding to the peak value. 

(pb) = 0.62 x 10 10 M q MpcT 3 at z — 0 (Planck Collabora- alisations, which is Gaussian when considering the logarithm 
tion et al. 2015). of the quantities. We estimated the parameter uncertainties as¬ 

suming that A and B are fully independent. This assumption is 
The parameter uncertainties in Eqs. 1 and 2 reflect the scat- probably not fully valid and thus the parameter uncertainties are 
ter among the power-law fits to the different Monte Carlo re- 
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PwHIM [10 Msun/Mp C ] PwHIM [^0 M sun /MpC ] 

Fig. 3. The WHIM density distribution at ^ld ~ 65 (where the WHIM density peaks at db ~26) using all the simulation data within the filamentary 
environments are shown with blue crosses and dashed lines. Red crosses indicate the distribution when using only volume elements with galaxies 
fainter than M r = -20.5 (left panel) or with gas hotter than T = 5 x 10 6 K 


only approximate. While the normalisation parameter A is for¬ 
mally constrained within 20%, the correlation between A and B 
(the higher normalisation may be compensated by a the steeper 
slope to some extent) is expected to yield larger variations for 
the WHIM density at given LD. To verify this, we calculated the 
68% confidence intervals for the WHIM density at different LD 
values using the WHIM density predictions using the power-law 
model with a pair of A and B values from each power-law fit to 
the randomised data (i.e. including the parameter correlations). 
In the range d h ~ 10-100 the WHIM density is uncertain by a fac¬ 
tor of 2-3, while at the lower densities (<5b —1—10) the variation 
reaches a factor of ~10 (see Fig. 5). 


4. Testing the method 

In the current work our emphasis is on a detailed description of 
our WHIM method. In a follow-up paper we will carry out a full 
testing of the reliability of our method using all available data 
of sufficient quality. Here we perform a first test on the accuracy 
of our LD-based WHIM column density predictions by compar¬ 
ing them with those obtained independently via X-ray absorption 
measurements of a background blazar. 

We chose to apply the first test to the Sculptor Wall (SW) and 
Pisces-Cetus (PC) superclusters, in the sight-line to the back¬ 
ground blazar H2356-309, since they are 1) among the closest 
reported X-ray absorption systems (z~0.03 and ~0.06) respec¬ 
tively (F10 and Z10) 2) covered with (2dF) galaxy survey data of 
sufficient quality (2dF, see Colless et al. 2001; Tago et al. 2006) 
whose 3) reported WHIM column densities are among the high¬ 
est in the literature. 

4. 1. Luminosity density fields and galactic filaments 

We applied the LD method (see Section 2.2) to the bj band lumi¬ 
nosities of the galaxies around SW and PC (see Fig. 6) and thus 
produced the LD/, fields (see Fig. 7). Around the redshifts of the 
X-ray absorption line centroids of both SW and PC, the blazar 
sight-line passes through enhanced LD regions (see Fig. 7). The 
enhancement of LD at the locations of the X-ray detected WHIM 
supports our finding from the simulations (Cl 2) that LD traces 
the WHIM. 


We then applied the Bisous model (see Section 2.1) to 
the above galaxy distribution to identify galactic filaments. We 
found that there are several filaments within the enhanced LD 
regions (see Fig. 6). There are 247 galaxies affecting the LD pro¬ 
files of SW and PC, i.e. located within two times the size of the 
LD field smoothing kernel (2.8 Mpc) from the H2356-309 sight¬ 
line at the matching redshifts. Only two of these were classified 
as non-filament galaxies. This supports our finding from simula¬ 
tions (Cl 2) that within the filamentary environments the WHIM 
density is enhanced from the cosmic average. Thus, it is justified 
to apply below the LD-WHIM density relation (Eq. 2) derived 
within the filamentary environments. 

Since both absorbers have originally been targeted due to a 
priori knowledge of the large-scale concentrations of galaxies 
around them, the above agreement serves as a positive indica¬ 
tor of the performance of our filament-finder. We will perform a 
more robust test in a follow-up work. 

The fortunate combination of the blazar position and the ori¬ 
entation of the filaments in both SW and PC results into rela¬ 
tively long (~10 Mpc) l.o.s. projection of high LD regions. This 
contributes to these systems being among the most significantly 
X-ray detected WHIM structures to date. 

4.2. Luminosity density profiles and N H 

In order to evaluate the column densities along the sight-line, we 
proceeded into characterising the radial behaviour of the lumi¬ 
nosity density. We performed this by producing the LD/, profiles 
of SW and PC by sampling the LD/, fields along the line of sight 
towards the blazar H2356-309 (see Fig. 8). By examining the LD 
fields and the LD profiles (see Figs. 7 and 8), we set a lower LD 
threshold LD/, jnln = 0.05 x 10 l(l L G Mpc for considering the 
radial profile of the LD field across the X-ray line centroid and 
set the likely redshift range of the absorber equal to that 

in which LD/,. > LD fcj , mi n . 

Within the reported uncertainty range of the redshift of 
the X-ray absorber (F10; Z10), there is a continuous, two- 
peaked structure of ~10 Mpc length in both SW and PC (see 
Fig. 8). The redshifts of the two peaks in the LD profiles are 
very similar. Their separation, A. ~ 0.001, corresponding to 
Aa ~ 0.01 A, at 20 A, is below the energy resolution of both 
Chandra/LETGS and XMM-Newton/RGS. Therefore we con- 
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Fig. 4. The distributions (blue crosses) of the power-law p whim = A x LDf normalisation A (left panel) and the index B (right panel) when fitting the 
randomised LD, -WHIM density relations. The best-fit log-normal models are shown with red solid line. The dotted red line shows the extrapolation 
of the best-fit model. The centroid and 68% intervals are indicated with the dashed and dotted vertical lines. 
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Fig. 5. The power-law approximation of the relation between the luminosity density (LD r ) and the WHIM density (p whim ) within the galaxy 
filaments in the simulations of C12 (solid line), and the l<x uncertainty interval (blue region). Left panel shows the full fitted LD range while the 
right panel shows the (S b —1—10 range. 


sidered the full ~10 Mpc long two-peaked structure as a single 
system for both SW and PC. 

We then converted our measured bj band Sculptor LD 
profiles into units of luminosity over-densities <5ld(z) = 

LD& (z)/(LD& ) (see the discussion in Section 3.2) by applying 
the mean bj band luminosity density in the full 2dF survey (see 
Table 1). We then used our LD-WHIM density relation (Eq. 2) 
to convert the above luminosity overdensity profiles into WHIM 
density profiles. Integration of the density profile then yielded 
the column density estimates using 

Afa(WHIM)= f 2 Pwhim(z) dz. (3) 

Jzi 

4.3. Uncertainties of the estimated WHIM column densities 

In this Section we analyse the two main sources of uncertainties 
affecting the estimates of the WHIM column densities. First, we 
analyse the reliability of the observed luminosity density pro¬ 


files. Second, we analyse the effect of the scatter in the calibra¬ 
tion of the LD-WHIM density relation in simulations (Cl 2). 


4.3.1. Uncertainties of the observed luminosity density field 

The main source of uncertainty in an observed LD profile is the 
shot noise due to the limited number of galaxies used for the 
LD field construction. We estimated this effect by applying a 
smoothed bootstrap procedure on the 2dF galaxy sample (see 
Liivamagi et al. 2012 for more details and justification of the use 
of the smoothed bootstrap in this context). The galaxies were 
selected with replacement as usual in bootstrap, but additionally 
their positions were randomised. The distance of the randomised 
position from the original one is defined by a smoothing kernel, 
whose size is half of that used for the original LD field smooth¬ 
ing in Section 2.2. Repeating this procedure 10000 times we 
obtained a set of bootstrapped LD profiles for SW and PC and 
we used the 68% LD scatter interval at each radius to determine 
the uncertainties of the LD profile measurements (see Fig. 8). 
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Fig. 6. The spatial distribution of 2dF galaxies around the H2356-309 sight-line (yellow line) projected at two orthogonal directions (upper and 
lower panels), at the distance range covering Sculptor Wall and Pisces-Cetus structures. The coordinates refer to the CMB rest frame. The galaxies 
belonging to filamentary environments are denoted with plus signs while other galaxies are denoted with crosses. The filament spines are indicated 
with green lines. The relative LD level at each radii is denoted with the blue line. The red parts of the sight-line highlight the radial ranges of the 
luminosity density profiles analysed in this work (see Fig. 8). 
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Fig. 7. The 2dF bj band luminosity density field slices of 1.4 Mpc thickness in two orthogonal directions (upper and lower panels) along the 
line-of-sight towards the blazar H2356-309 (red horizontal lines) in the Sculptor Wall (left panels) and Pisces-Cetus (right panels) structures. The 
contours indicate our adopted lower LD threshold LD/, ;>m i n = 0.05 x 10 10 L o Mpc~ 3 for extracting the data for further analysis. The radial ranges 
equal those used for the luminosity density profile analysis (see Fig. 8). The galaxies relevant to the sight-line, i.e. within a distance of two times the 
smoothing kernel size from the sight-line (2.8 Mpc), are marked with dark red symbols: plus signs correspond to locations of galaxies belonging 
to a filamentary environment while other galaxies are denoted with crosses. The purple vertical lines indicate the centroids of the Chandra X-ray 
absorption lines (F10; Z10). 


We propagated this scatter into N# measurement by convert¬ 
ing each bootstrapped LD profile into a WHIM density profile 
using Eq. 2 and integrating the profiles using Eq. 3. The resulting 
N h distributions followed well the log-normal prediction (see 
Fig. 9). We thus used the best-fit log-normal centroid and width 
parameters to obtain the best values and lcr intervals of N h i.e. 


0.8 [0.6-1.1] x 10 20 cm- 2 and 0.5 [0.4—0.7] x 10 20 cm- 2 for SW 
and PC, respectively. 

4.3.2. Effect of scatter in the LD-WHIM density relation 

Due to the significant scatter of the WHIM density values for 
a given LD value in the cosmological simulations used in this 
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Redshift in heliocentric frame 
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Fig. 8. The 2dF bj band luminosity density profiles of the Sculptor Wall (left panel) and Pisces-Cetus (right panel) along the H2356-309 sight-line 
(black solid lines) together with lcr uncertainties (blue regions). The horizontal (green) dashed line indicates our adopted lower LD threshold 
LD= 0.05 x 10 10 L 0 Mpc~ 3 for extracting the data for further analysis. The vertical dashed lines indicate the best-fit heliocentric redshifts 
of the Chandra X-ray lines while the dotted lines bracket the statistical 1 cr uncertainties of the redshifts (from F10 and Z10). The co-moving 
distances refer to the CMB rest frame. 



19.4 19.6 19.8 20.0 20.2 

log N h (cm -2 ) 



log N h (cm -2 ) 


Fig. 9. The WHIM hydrogen column density distributions (black crosses) due to the LD profile measurement uncertainties for the Sculptor Wall 
(left panel) and Pisces-Cetus (right panel) structures. The best-fit log-normal models are shown with (red) solid lines. The best value and the lcr 
confidence intervals are indicated with (green) dashed line and the (blue) shaded regions, respectively. 


work, our model prediction is statistical by nature. The following 
list contains brief discussion about the possible origins of the 
scatter. 

- The galaxy luminosity and the WHIM density do not solely 
depend on the underlying dark matter density but may de¬ 
pend on e.g. the stellar feedback mechanisms and the gas 
physics. 

- The HOD used to populate the dark matter halos with galax¬ 
ies is statistical by nature. 


- The dark matter halo geometry may deviate from the as¬ 
sumed spherical symmetry in the HOD modelling. 

We already characterised this scatter in Section 3.5 by the 
power-law fits to the set of randomised LD-WHIM density rela¬ 
tions, when deriving the uncertainties of the LD-WHIM density 
power-law model parameters. Due to the parameter correlations 
(see Section 3.5) we did not perform a simple analytical error 
propagation using Eqs. 2 and 3. Rather, we used here the pairs 
of the best-fit power-law normalisation and index parameters ob- 
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Fig. 10. The WHIM hydrogen column density distributions (black crosses) due to the LD-WHIM density scatter in the simulations (Cl 2) for the 
Sculptor Wall (left panel) and Pisces-Cetus (right panel) structures. The best-fit log-normal models (red lines) in the fitted range (solid) are shown 
together with the extrapolation (dotted). The best value and the lcr confidence intervals are indicated with (green) dashed line and the blue shaded 
regions, respectively. 


tained in Section 3.5 to obtain a set of LD-WHIM density rela¬ 
tions which reflects the scatter of the WHIM density values for a 
given LD in C12 simulations. We applied each of these scattered 
LD-WHIM density relations to the measured LD profiles of SW 
and PC first to obtain a set of WHIM density profiles and then, 
by integration, we obtained distributions of Nh values for both 
structures (see Fig. 10). 

At Nh values within the 68% confidence interval around the 
centroid, the distributions can be reasonably well approximated 
with a log-normal model. We thus used the best-fit log-normal 
centroid and width parameters to obtain the best values and lcr 
intervals of Nh i.e. 0.8 [0.6-1.1] x 10 20 cm~ 2 * and 0.5 [0.4— 
0.8] x 10 20 cm~ 2 for SW and PC, respectively. At the lowest and 
the highest Nh values the distribution significantly exceeds the 
log-normal distribution. The low end tail is a consequence of the 
high end tail of the index distribution (see Fig. 4). Note that the 
rather large (by a factor of ~10) uncertainty of the LD-WHIM 
density relation at the lowest densities (see Fig. 5) has a negligi¬ 
ble effect in case of SW and PC, since most of the measured LD 
values are much higher (see Fig. 8). 

4.4. Final N H values 

The effects of 1) the Sculptor LD profile measurement uncer¬ 
tainties using the 2dF data and 2) the LD-WHIM density scatter 
in the cosmological simulations, on the Nh estimates are similar 
and significant (i.e. a variation by a factor of ~2). We thus com¬ 
bined the two uncertainties in quadrature. The best values differ 
slightly and thus we averaged those, obtaining the final results 
that are given in Table 2. 

4.5. Comparison with X-rays 

Our procedure predicts the hydrogen column density of the 
WHIM from the galaxy luminosity density. However, the typical 
signatures of the WHIM in the X-ray band are line features from 
highly ionised metals, like O vn. The conversion between the hy¬ 


drogen column density and metal ion column density depends on 
the ionisation fraction (i.e. the temperature) and the metal abun¬ 
dance. If only a single absorption line is measured, as in the case 
of the Sculptor Wall structure (F10), the temperature and col¬ 
umn density constrains are very poor. If we assume that oxygen 
abundance is 0.1 Solar and that the temperature is 10 6 K (i.e. the 
ionisation fraction of Ovu is 1.0) our Nh estimate for the SW 
structure (log Ah = 19.9^'j) translates to log Aovn = 15.7;[q 2 , 
consistent within the large uncertainties of the X-ray measure¬ 
ment of FI 0, log Aovh = 16.8^qq (see Table 2). 

In the case of the Pisces-Cetus structure, lines from several 
elements and ionisation stages were reported (Z10) for a warm 
absorber 2 (log T (K) ~ 5.4). Thus the temperature and the 

equivalent hydrogen column density were well constrained and 
we can make a direct comparison (assuming that oxygen abun¬ 
dance is 0.1 Solar). Our estimate (log Ah = 19.7^j' 2 ) is con¬ 
sistent within lcr uncertainties with that of Z10, i.e. log Ah = 
20 . 1 + 0 . 2 . 

The agreement of our LD-based prediction for the column 
density of the WHIM in SW and PC structures with that obtained 
from X-ray absorption measurements indicates that our method 
is robust. This also indicates that the luminosity density and the 
galactic filaments are reliable signposts for the WHIM. 

5. Galaxy confusion 

In the above we have interpreted the X-ray absorption features 
being due to warm-hot intergalcictic medium. We examine here 
an alternative hypothesis, the galaxy confusion, which states 
that instead of the intergalactic diffuse WHIM, the absorption 
is caused by the X-ray halos of accidental galaxies close to the 
studied sight-line (e.g. Williams et al. 2013). Note that in this 
case the match of our intergalactic LD - based WHIM column 
densities and the X-ray ones (F10; Z10) is not causal, but rather 
a co-incidence. 

2 Z10 also reported a possible hotter component at the same redshifts 

but due to its lower significance we do not consider it here. 
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Table 2. Properties of the absorbers. 



X-ray 

LD 

X-ray 

LD 


Sculptor Wall 

Pisces-Cetus 

z° 

0.0310-0.0334 

0.0286-0.0320 

0.0618-0.0628 

0.0608-0.0630 

N h ( 10 2 °cnr 2 ) 

- 

0 . 8 + 0 .4 

- 

0.5+0.3 

log N„ 

- 

19 ' 9 -Q3 

20 . 1 ± 0 . 2 & 

19 7 + 0.2 

Ly ‘ -0.3 

log N 0 v ii c 

16.8^- 3 

15 7 +9 - 2 

_ iJ -o.3_ 

- 



Notes. (a) X-ray: The uncertainty interval of the Chandra X-ray absorption line redshift in Sculptor Wall (F10) and Pisces-Cetus (Z10), LD: The 
extent ( Z 1 -Z 2 ) of the LD structure intersected by the H2356-309 sight-line. The redshifts refer to the heliocentric frame. <b) The more significant 
warm component <c) Derived from assuming T = 10 6 K, O abundance of 0.1 Solar, and Solar O/H ratio from Grevesse & Sauval (1998). 


Table 3. Nearby galaxies 


ID 

RA 

DEC 

z“ 


M bj d b type c 


deg 

deg 


O 

© 

kpc 


Sculptor Wall 


1 

359.79297 

-30.58761 

0.0296 

0.08 

-17.3 

90 

star forming 

2 

359.70248 

-30.56833 

0.0297 

0.04 

-16.7 

190 

star forming 

3 

359.90756 

-30.65583 

0.0301 

0.62 

-19.5 

240 

elliptical 

4 

359.84072 

-30.80313 

0.0295 

0.29 

-18.7 

390 

low star formation 

5 

0.12302 

-30.72346 

0.0316 

0.08 

-17.3 

700 

star forming 

6 

359.30189 

-30.61506 

0.0295 

0.02 

-15.8 

880 

low star formation 

7 

359.36436 

-30.46063 

0.0312 

0.88 

-19.9 

890 

low star formation 

8 

0.27432 

-30.61173 

0.0316 

0.06 

-16.9 

970 

star forming 


Pisces-Cetus 


9 

359.79467 

-30.62261 

0.0627 

0.14 

-17.9 

50 

active starburst 

10 

359.87066 

-30.55599 

0.0625 

0.21 

-18.4 

470 

active starburst 

11 

359.79755 

-30.76000 

0.0613 

0.13 

-17.8 

590 

star forming 

12 

359.60680 

-30.64623 

0.0612 

0.09 

-17.4 

680 

low star formation 

13 

359.76432 

-30.84869 

0.0623 

0.11 

-17.6 

990 

low star formation 


Notes. un Redshifts are reported in the heliocentric frame. (i,) distance from the H2356-309 sight-line in the plane of the sky. (cl Based on Madgwick 
et al. (2002). 


The sole existence of a galaxy within a virial radius from 
an X-ray absorber is obviously no proof that the galaxy and 
its circum-galactic gas are the origin of the X-ray absorption 
features that we regard as WHIM signatures. We explored the 
galaxy confusion hypothesis quantitatively in this work by 1 ) es¬ 
timating the density of the X-ray halos of the nearby galaxies at 
the distance of H2356-309 sight-line and 2) comparing these es¬ 
timates with the X-ray measurements. As possible candidates we 
very generously considered all such 2dF galaxies within the red- 
shift range matching those of the SW and PC structures, which 
are located within 1 Mpc (i.e. several times the virial radii) from 
the H2356-309 sight-line. There are eight (five) of such candi¬ 
dates for SW (PC), (see Table 3). 

In the above list there are no massive galaxies which could 
be obvious candidates for the galaxy confusion. It is possible 
that close to the sight-line there are galaxies fainter than the 2dF 

limit, corresponding to Mb -16 at these redshifts. However, 

these luminosities correspond to ~1% of that of the Milky Way. 
Thus, one needs 100 of such objects to reach the mass of the 
Milky Way and possibly a similar WHIM K H level of 10 19 cm -3 


(assuming the temperatures of such haloes do reach WHIM val¬ 
ues). These objects should be clustered at the sight-line, at the 
redshifts consistent with the X-ray measurements. Thus these 
objects should be more clustered than Milky Way - like galaxies 
(see Table 3), i.e. their bias parameter should be larger. This is in¬ 
consistent with observations and simulations that show that faint 
galaxies are less clustered than the luminous ones (e.g. Norberg 
et al. ( 2002 a); Cl 2 ). 

Galaxy 3 in our list is the only elliptical galaxy while all oth¬ 
ers are spirals. In the following we examine separately the two 
galaxy types. For the comparison with other works on galaxies, 
we approximated the absolute bj band magnitudes and luminosi¬ 
ties with those in the /(-band (Norberg et al. 2002b). 

5. 1. Elliptical galaxy 3 

Williams et al. (2013) suggested that the elliptical galaxy 3 in 
our list, at a distance of 240 kpc from the H2356-309 sight¬ 
line at the redshift of the Sculptor Wall, may be responsible 
for the X-ray absorption measured by F10. The bj band mag- 
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nitude of the galaxy corresponds to Lg ~ 10 10 Lg 0 - From the 
analysis of X-ray halos of elliptical galaxies (Fukazawa et al. 
2006), the corresponding X-ray luminosity of this object would 
be Lx ~ 10 4() erg s _1 and a temperature of T~ 10 6 K, consis¬ 
tent with the WHIM properties. Extrapolating the correspond¬ 
ing gas density profile in (Fukazawa et al. 2006) to a radius of 
240 kpc (corresponding to the impact parameter of the sight-line 
to H2356-309) yields a density level of ~ 10 ~ 5 cm -3 . With such 
a low density a path length of ~ 10 25 cm, i.e. ~ 10 Mpc would 
be required to build up the measured level of the column den¬ 
sity, i.e. N h ~ 10 20 cm -2 . This is clearly inconsistent with the 
virial radius of the galaxy, estimated as 350 kpc in Williams et al. 
(2013). 

5.2. Spirals 

All the spiral galaxies in the candidate list have a much smaller 
or similar B-band luminosity Lg compared to the Milky Way, 

whose M b ~ M*.-20 (Binney & Merrifield 1998; Pritchet & 

van den Bergh 1999; van den Bergh 2000). We wish to use the L B 
values for comparing the stellar masses of the candidate spirals 
with that of the Milky Way. This is complicated since the /1-band 
may be “contaminated” due to star formation. All our candidate 
spirals are star-forming and thus they likely have a similar bias 
in the Z?-band - based mass estimates of the candidates and the 
Milky Way. However, in order to form a conservative estimate 
we assumed a factor of 100 scatter in the mass /L B ratio, render¬ 
ing the most massive candidate spirals 100 times more massive 
than the Milky Way. 

The intensity of the X-ray emission from T ~ 10 6 K gas in 
several edge-on spiral galaxies, e.g. NGC 5775 (Li et al. 2008) 
and NGC 4631 (Yamasaki et al. 2009), has been observed to 
decrease exponentially with the height from the galactic plane. 
The joint X-ray emission and absorption measurements of O vn 
and O vm lines of the Milky Way halo have obtained tempera¬ 
tures ~10 6 K (consistent with WHIM) and equivalent hydrogen 
column densities at 10 19 cm ~ 2 level (e.g. Hagihara et al. 2010; 
Sakai et al. 2014), i.e. ~ 10% of that consistently determined by 
us and via X-rays in Sculptor Wall and Pisces-Cetus. The result¬ 
ing path length scale of the absorber and thus the extent of the 
WHIM halo in spiral galaxies in the above works is consistently 
constrained below 10 kpc. While this indicates that the halo is 
identified with the stellar component, this also conflicts with the 
reports of the Galactic ~100 kpc WHIM halo (e.g. Gupta et al. 
2012). In fact, Wang & Yao (2012) reported numerous signifi¬ 
cant errors in Gupta et al. (2012) which resulted in an order of 
magnitude overestimation of the Galactic WHIM halo size. 

Since our candidate spirals have masses less than 100 times 
that of the Milky Way (see above), we assume here conserva¬ 
tively that all our candidate spirals have WHIM halos with a 
similar exponentially decreasing density profile as that of the 
Milky Way, except that scaled up by a factor of 100, i.e. with 
a central density of 10 1 cm ~ 2 and a scale height of 3 kpc. The 
most likely candidate for galaxy confusion is galaxy 9 (see Ta¬ 
ble 3), the nearest one to the absorbers. Williams et al. (2013) 
associated this galaxy with the X-ray absorber at Pisces-Cetus 
(Z10). Extrapolating the exponential density profile to 50 kpc 
(the distance between the galaxy and H2356-309 sight-line), we 
obtain a very low WHIM density of 10 8 cm -3 . Let’s assume 
that the column density of ~ 10 20 cm ~ 2 is achieved along an ab¬ 
sorber path length of 10 kpc (~ 10 22 cm). Consequently, the re¬ 
quired number density of the WHIM halo is ~ 10 2 cm ~ 3 at the 
distance of the H2356-309 sight-line from the centre of a given 
galaxy. Thus, the exponential density profile of galaxy 9 fails 


totally to produce the measured level of WHIM column density. 
Since the other galaxies are more distant, the discrepancy is even 
bigger for those. Unless we underestimated the candidate spiral 
to Milky Way mass ratio by a factor of 10 6 due to the star for¬ 
mation issue, the candidate spirals cannot produce the measured 
absorption in Pisces-Cetus. 

In general, considering the projection effects, a Milky Way - 
like galaxy must be located within ~ 10 kpc to the blazar sight¬ 
line in order to produce WHIM column densities at a level of 
~ 10 19 cm -2 . 


6. Discussion 

In this work we tested the ability of our method to infer the pres¬ 
ence of the WHIM by targeting large scale structures associated 
to putative WHIM detections. In the future we shall apply our 
method to perform a blind WHIM search along the sight-lines 
to distant quasars located within the areas of the 2dF and SDSS 
galaxy surveys. Note that our method is not limited to the direc¬ 
tions towards the brightest background blazars, as is the case for 
the current rare WHIM detections based on the identification of 
characteristic X-ray and far ultra violet absorption lines. Rather, 
our method has potential for extending the hunt for the missing 
baryons to the full volumes of the spectroscopic galaxy surveys 
like SDSS and 2dF. 

Our method can be used to estimate the mass density of 
the WHIM. Combined with the independent X-ray or FUV ab¬ 
sorption measurements of the same WHIM structures, our ad¬ 
ditional constrain has potential for improving the diagnostics of 
the WHIM physics. In particular, there is a possibility of break¬ 
ing the degeneracy between the WHIM column density and the 
abundance of the absorbing metal, and thus providing constraints 
for the latter. This could improve our understanding of the chem¬ 
ical and thermal enrichment of the intergalactic medium via su¬ 
pernovae explosions and AGN feedback. Currently these effects 
are poorly known and modelled in the hydrodynamical simula¬ 
tions. 

Our LD-WHIM density relation together with the filament¬ 
tracing technique thus has potential for mapping a large sam¬ 
ple of WHIM properties with the currently existing data. This 
is necessary for properly assessing the contribution of WHIM 
to the cosmological mass density budget. Thus our method has 
potential for advancement in the cosmological missing baryons 
problem in the near future. Our methods and the potential WHIM 
maps could be used to optimise the observational strategies for 
the search of the missing baryons with next generation satellites 
like Athena. 


7. Conclusions 

In this work we proposed, implemented, applied and tested a 
novel method to identify the elusive WHIM using galaxies in 
the filaments of the cosmic web. Our method rests upon a corre¬ 
lation between the density of the “invisible” WHIM and that of 
the observed galaxy luminosity. The calibration of the method 
was performed using the hydrodynamical numerical simulation 
(Cl2) and HOD-based mock galaxy catalogues mimicking the 
properties of SDSS galaxies, following the Zehavi et al. (2011) 
prescriptions. To test the performance of our method on observa¬ 
tional data we applied it to the distribution of galaxies in the 2dF 
survey, focusing on the Sculptor Wall and Pisces-Cetus super¬ 
clusters. The main results, obtained both in the calibration and 
application phases are summarised below. 
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- In the C12 simulation the WHIM is preferentially associated 
to the filamentary structures. We identified the latter by ap¬ 
plying the Bisous model (Stoica et al. 2005; Tempel et al. 
2014a) to the distribution of the mock galaxies and found 
that the mass fraction of the gas in WHIM phase associated 
to these structures is ~70%, by a factor of ~1.6 higher than 
the average in the full simulation box. 

- When considering the WHIM gas and the galaxy lumi¬ 
nosity in the filament environments, we found a signifi¬ 
cant correlation between their densities. The Pearson corre¬ 
lation coefficient between the two turned out to be ~80% 
and the relation between the WHIM gas overdensity and 
the galaxy luminosity overdensity consistent with linear: 
<5whim = 0.7 + 0.1 x 5 ld°' 9±0 ' 2 - This suggests that both 
the diffuse gas in the WHIM and the stellar component in 
the galaxies trace the same underlying dark matter density 
field. 

- The above relation can be used in reverse, i.e. given the ob¬ 
served luminosity density in the filamentary structures iden¬ 
tified in the spatial distribution of galaxies in galaxy redshift 
surveys one can use the relation to infer the column density 
of the WHIM. 

- When testing our method on observational data, we found 
that the luminosity density was significantly enhanced from 
the average at the locations of independently detected WHIM 
(F10; Z10). This indicates that the luminosity density traces 
the WHIM, consistently with what we found in the simula¬ 
tions. 

- Our LD-based predictions for the column density of the 
WHIM in the Sculptor Wall and Pisces-Cetus superclusters 
agree with those obtained via X-ray absorption. This agree¬ 
ment indicates that our method is robust in estimating the 
density of the WHIM. Also, the galaxy filaments and the lu¬ 
minosity density are reliable signposts of the WHIM. 

- The fortunate combination of the angular position of the 
bright background blazar H2356-309 behind suitably ori¬ 
ented, long filaments in the Sculptor Wall and Pisces-Cetus 
superclusters results into relatively long (~ 10 Mpc) line-of- 
sight projection of high LD regions. This contributes to these 
systems being among the most significantly X-ray detected 
WHIM structures to date. 

- The signal that we detected cannot originate from the halos 
of the nearby galaxies since they cannot account for the large 
WHIM column densities that our method and X-ray analysis 
consistently find in the Sculptor Wall and Pisces-Cetus su¬ 
perclusters. 
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Fig. A.l. The halo occupation number, i.e. the mean number of galax¬ 
ies in a halo of a given virial mass, is shown in different luminosity bins 
(colour coding as in Fig. A.3). The solid histograms represent the statis¬ 
tical realisations of the galaxies that populate the halos extracted from 
the (Cl 2) simulation and their substructures. The dotted curve shows the 
predictions from Eq. A.2 using the parameters in Zehavi et al. (2011). 


Appendix A: HOD Formalism 

We applied halo occupation distribution (HOD) formalism to our 
adopted large-scale simulation of Cui et al. (2012) in order to 
construct the galaxy distribution. Dark matter halos and subha¬ 
los represent tracers of the central galaxies and satellites, respec¬ 
tively (Peacock & Smith 2000; Kravtsov et al. 2004; van den 
Bosch et al. 2007). We thus started by producing a DM halo cat¬ 
alogue by applying the Friend-of-Friend (FoF) algorithm with 
linking length parameter b = 0.2 (in the units of the mean inter¬ 
particle separation) to the simulated data. This choice of the link¬ 
ing length value yields mass functions consistent with different 
theoretical predictions of the virial mass function (e.g. Gao et al. 
2004; Springel 2005). The resulting FoF catalogue contains the 
position and the virial mass of each DM halo. In the following 
we describe how we used the FoF catalogue to obtain the mag¬ 
nitudes and positions of the galaxies. 

1. For each DM halo we used its FoF mass from above to cal¬ 
culate the virial radius R vlr . to know up to which scale to 
populate it with subhalos. This was done according to the 
spherical collapse formalism which yields the estimate of the 
virial overdensity A vlr (z) (e.g. Peebles 1980; Eke et al. 1996; 
Kitayama & Suto 1996; Bryan & Norman 1998) which is 
linked to the virial mass and the virial radius as: 

Mh = (A.l) 

3 12 m (z) 

where p c (z ) and Q m (z) represent the critical density and the 
cosmological matter density parameter at redshift z, respec¬ 
tively. 

2. We assigned each halo with a concentration parameter ac¬ 
cording to mass-concentration relation of Bullock et al. 
(2001), assuming a log-normal scatter cr lne = 0.25 around 
the mean value. 



subhalos (dotted line and dots) hosting galaxies with different r band 
absolute magnitudes, when applying the HOD formalism to C12 simu¬ 
lations (colour coding as in Fig. A.3). 


3. We populated each halo with subhalos by performing Monte- 
Carlo realisations of the subhalo mass function model of 
Giocoli et al. (2010) which features both a redshift evolu¬ 
tion and a concentration dependence of the subhalo mass 
function. The model assumes that the spatial distribution of 
the subhalos is less concentrated than the NFW DM profile 
(Navarro et al. 1996), since this includes the effects of the 
dynamical friction and the tidal stripping (Gao et al. 2004; 
van den Bosch et al. 2004; Giocoli et al. 2008). Thus, for 
each halo we now have the subhalo populations with known 
positions and masses. 

4. We then assigned each halo and subhalo with a galaxy with 
a luminosity value according to the abundance matching ap¬ 
proach (see e.g. Behroozi et al. 2010) as follows. In the core 
of the method is the halo occupation function 


(N(M h )) = - 


! / log M/, - log M„ 

\ trlog m 


X 


1 + 


i M h -M 0 y 


(A.2) 


which describes the mean number of halos within a parent 
halo of mass M/, (Zehavi et al. 2011). 

We then made the standard assumptions that 1) all subha¬ 
los are populated by one galaxy only, i.e. that the number of 
substructures into which we have resolved the parent halos is 
sufficient to host at most one galaxy and that 2) there are no 
“orphan” subhalos. The latter hypothesis is justified by the 
fact that there is no evidence for massive dark halos with no 
baryon content. Using these assumptions Eq. A.2 then de¬ 
scribes the mean number of galaxies within a parent halo of 
mass Mi,. 

The values of the parameters of Eq. A.2 have been found by 
Zehavi et al. (2011) by fitting the projected SDSS-DR7 2- 
point galaxy-galaxy correlation function of Abazajian et al. 
(2009) in the luminosity range M r = -[18.5,-22] sam¬ 
pled in seven, equally spaced bins. Thus, the parameters of 
Eq. A.2 are different for the different luminosity bins. 
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We then applied the method to our FoF DM halos obtained 
from the C12 simulations as follows: for a given parent halo 
of mass M/, we use Eq. A.2 to determine the mean number of 
galaxies at a given magnitude M r (see Figs. A. 1 and A.2). We 
repeated the procedure for each magnitude bin from M r = 
-18.5 to M r = -22 thus obtaining the luminosity function 
N (L r ) for a given host halo. 

We then ordered the above galaxies from most luminous to 
least luminous. For the same parent halo we went back to 
the subhalo distribution obtained above (Giocoli et al. 2010) 
and ordered the subhalos from the most massive to the least 
massive. We then matched the parent halo or its subhalo of 
a given mass ranking by a galaxy with the same luminosity 
ranking (most massive with the most luminous etc.) until all 
galaxies were assigned to the subhalos. We removed the low 
mass subhalos which were assigned with no galaxies. 

The outcome is a set of parent halos, extracted from the 
C21 simulation, containing a set of subhalos. Each halo and 
subhalo has a galaxy at its centre with known location and 
r-band luminosity (see Fig. A. 3). 
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Fig. A.3. Three orthogonal projections of the distribution of satellite galaxies in centres of DM halos of ~1 Mpc radius in our adopted simulation 
of C12 (coloured dots). The colour coding indicates the magnitude of a given galaxy. The black dots show the positions of galaxies fainter than 
M,.=-18 that populate halo and sub-halos in the simulation according to the mass function by Giocoli et al. (2010) down to 10 10 M e . 
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